Uniaxial Strain in Graphene by Raman Spectroscopy: G peak splitting, Gruneisen 

Parameters and Sample Orientation 



T. M. G. Mohiuddin^ A. Lombardo^ R. R. Nair\ A. Bonetti^ G. Savini^ R. Jalil^ 
N. Bonini^, D.M. Basko C. Galiotis^, N. Marzari^, K. S. Novoselov^ 
A. K. Geimi, A. C. Ferrari^l 

^Department of Physics and Astronomy, Manchester University, Manchester UK 
^Engineering Department, Cambridge University, Cambridge, UK 
Department of Materials Science and Engineering, 
Massachusetts Institute of Technology, Cambridge, MA, USA 
Laboratoire de Physique et Modelisation des Mileux Condenses, 
Universite Joseph Fourier and CNRS, Grenoble, France 
^ FORTH/ ICEHT and Materials Science Dept., University of Patras, Patras, Greece 

Graphene is the two-dimensional building block for carbon allotropes of every other dimensionality. 
Since its experimental discovery, graphene continues to attract enormous interest, in particular as 
a new kind of matter, in which electron transport is governed by a Dirac-like wave equation, and 
as a model system for studying electronic and phonon properties of other, more complex, graphitic 
materialsP,[3,i,ll. Here, we uncover the constitutive relation of graphene and probe new physics of 
its optical phonons, by studying its Raman spectrum as a function of uniaxial strain. We find that 
the doubly degenerate E29 optical mode splits in two components, one polarized along the strain and 
the other perpendicular to it. This leads to the splitting of the G peak into two bands, which we call 
G"*" and G~, by analogy with the effect of curvature on the nanotube G peak[l,0,|3|. Both peaks red 
shift with increasing strain, and their splitting increases, in excellent agreement with first-principles 
calculations. Their relative intensities are found to depend on light polarization, which provides 
a useful tool to probe the graphene crystallographic orientation with respect to the strain. The 
singly degenerate 2D and 2D' bands also red shift, but do not split for small strains. We study the 
Gruneisen parameters for the phonons responsible for the G, D and D' peaks. These can be used to 
measure the amount of uniaxial or biaxial strain, providing a fundamental tool for nanoelectronics, 
where strain monitoring is of paramount importance[l, @| 



Strain arises when a crystal is compressed or stretched 
out of its equilibrium shape, with the stiffness ten- 
sor providing the constitutive relation between applied 
stress and final strain state. Atomic relaxations of- 
ten accompany the process, also resulting in an effec- 
tive renormalization of the constitutive relations. The 
presence of strain can significantly affect device perfor- 
mance. Sometimes, strain is intentionally applied to 
improve mobility, as in the strained silicon technology, 
which is used in modern microelectronics. Thus, the pre- 
cise determination and monitoring of stress and strain 
is a key requirement [1, 0|- Strain modifies the crystal 
phonons, with tensile strain usually resulting in mode 
softening, and the opposite for compressive. The rate 
of these changes is summarized in the Gruneisen pa- 
rameters, which also determine the thermomechanical 
properties Thus, monitoring phonons is often the 
clearest and simplest way to detect strain and, if the 
Gruneisen parameters are known, to quantify it. 

Raman spectroscopy has emerged as the main tech- 
nique to probe graphene's phonons[ll|. It can iden- 
tify the number of layers in a samplefllj. determine the 
amount of doping and presence of disorder 
[l6| . study graphene's edges 17, 13, 3, 2ol . 

tanharmonic processes and thermal conductivity [2^ 
Raman studies of graphene also revealed novel phys- 
ical phenomena, such as Kohn anomalies [H], and the 



breakdown of the Born-Oppenheimer approximation [12, 
Idjiil, 0|- In all these cases, experimental observations 



2l| and quan- 



have successfully partnered with first-principles calcu- 
lations, the latter providing additional microscopic in- 
sights and understanding, while being validated by the 
comparison with measurements. The Gruneisen param- 
eters for the vibrational modes of graphite under bi- 
axial strain were calculated by first-principles, yielding 
excellent agreement with the thermomechanical prop- 
erties of graphitefi^l- Recently, changes to the Ra- 
man spectra were reported due to the presence of stress 
in graphenefil S S, llil, [H, but the inferred 



strains disagreed by a factor 5 or more for similar Ra- 



man shifts^, \m 31 



l32l |. Furthermore, no significant 
difference was seen between the cases of uniaxial and bi- 
axial strain 2^ 31 , S^l , in contrast with theory, and the 
opening of a band gap at the K point was suggested [Hi, 
again in contrast with theory for small strains. It is thus 
necessary to conduct an accurate study in order to un- 
cover the physics of strain for the graphene phonons. 

In this work, we carefully apply uniaxial strain up to 
^--^1.3% to a graphene monolayer, in typical steps of 0.05% 
(minimum step 0.01%; maximum 0.25%) using two and 
four point bending setups as described in Methods (see 
Fig. 1), and compare this with first-principles calcula- 
tions. The Raman spectra measured at each step are 
fully reproducible over multiple loading and unloading 
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Figure 2: (Color Online) (a) G (b) 2D peaks as a function of 
uniaxial strain. The spectra are measured with incident light 
polarized along the strain direction, collecting the scattered 
light with no analyzer, see Methods. Note that the doubly 
degenerate G peak splits in two subbands, G+ and G~, while 
this does not happen for the 2D peak. The strains, ranging 
from to ~0.8%, are indicated on the right side of the spectra 



Figure 1: (Color Online) (A) Scheme (not to scale) of the 
substrate coated with SU8. A graphene monolayer is placed 
in the middle; (B,C) Scheme (not to scale) of (B) two point, 
and (C) four point bending set up. Note that a typical sample 
is 10'*-10'* smaller than the substrate length, see Methods. 



cycles, with no hysteresis. This allows us to clarify the 
picture for Raman spectra in strained graphene. 

Figure 2 plots some representative spectra as a func- 
tion of strain. The origin of the main Raman peaks is ex- 
plained in Methods. The strain is parallel to the longest 
side of the substrate (Fig 1), and is given by the ratio 
of substrate thickness to twice the radius of curvature. 
The spectra are fitted with lorentzians, and Fig. 3 plots 
the resulting trends for the G and 2D peaks. Note that 
Figs. 3a, b are a combination of over 80 measurements on 
two samples, strained in two different experimental set- 
ups, and include a loading, unloading and final loading 
cycle. Within the spectrometer resolution we find no dif- 
ference on pre-history and, for a single sample and cycle, 
the strain dependence is smooth. Linear fits using all the 
data yield duc+lde —10.8 cm'V %; duc-lde —31.7 
cm^V %; duj2D/ds —64 cm'^ / % and du2D'/de —35 
cm~^/ %. Where we call G+ and G~ the higher and 
lower G sub-bands, by analogy with nanotubesjl, 0|- 

The observed behavior can be explained by considering 



the effect of uniaxial strain on the optical modes respon- 
sible for the G, and D and D' peaks, respectively. 

The Griineisen parameter for the doubly-degenerate, 
in-plane, Raman active i?2g phonon, 'yE2g, is[10l|: 



IE2 



den 



(1) 



where en = £11 + £tt is the hydrostatic component of the 
applied uniaxial strain, I is the longitudinal direction, 
parallel to the strain, and t is the direction transverse to 
it; uj%^ is the G peak position at zero strain. The shear 



deformation potential /?_E2g is defined as|34l. l35 

PE2„ = — 



d£s 



(2) 



where Ss = £11 — £tt is the shear component of the strain. 

Under uniaxial strain, the solution of the secular equa- 
tion for the E2g mode isjsi, [H, \M, 



,7-B2g {£u+£tt) ± ^Pe2,(^E2 



{eu 



£tt) 
(3) 



where Aw^^^ is the shift resulting from the hydrostatic 
component of the strain, and A.uj% is the mode split- 
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Figure 3: (Color Online) Positions (Pos) of the (a) G"'" and 
G^, and (b) 2D peaks, as a function of applied uniaxial strain. 
The blue lines are linear fits to the data. The slopes of the 
fitting lines are also indicated 



nanotubes will have a TO G~ peak and a LO G"*", with 
the former softer than the latter, and more sensitive to 
diameter changes. This simple picture is reasonable for 
semiconducting nanotubes[5|, while in metaUic, a fur- 
ther significant softening of the LO mode takes place due 
to the enhanced Kohn anomaly resulting from electron 
confinement However, this further effect must be ab- 
sent in "unrolled" tubes, i.e. graphene. Indeed, the Full 
Width at Half Maximum (FWHM) of the G+ and G^ 
peaks in graphene is roughly constant as a function of 
strain at ~ 12cm~^, whereas FWHM(G~) in metallic 
nanotubes becomes much larger, due to the increased 
electron-phonon coupling contribution 

By fitting the trends in Fig. 3 to Eqs. 1,2, 3 we 
can experimentally determine the Gruneisen parameters 
for graphene. Under uniaxial strain[l3| £u = £ and 
Ett = —ve. Where u is the Poisson's ratio. If one could 
strain free-hanging graphene samples, the Poisson's ratio 
for graphene itself should be used. This can be taken as 
the in-plane Poisson's ratio of graphite '^0.13 [Hi. How- 
ever, the lack of loading-unloading hysteresis for our re- 
sults implies good adhesion between graphene and our 
substrates for the whole range of applied strains. SU8 
is a transversely isotropic material with a 0.33 in plane 
Poisson's ratio[39j. PET and perspex have also Poisson's 
ratios between 0.3-0.35. We thus use i^=0.33. This cor- 
responds to the case of ideal contact between graphene 
and substrate. Eq. 3 is now rewritten as: 



Au 



yielding: 
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2uJGo (1 - £ 



(4) 



(5) 



ting due to the shear component of the strain. All)q+ = 
^^E^g and AujQ- = Auj]^^^ are the shifts of the G+ and 
G~ peaks relative to zero strain. 

It is important to note that the resulti ng p honon eigen- 
vectors are orthogonal to each other sl, 33, 3^,[33|, with 



/3. 



Auj, 



E2„ 



G+ 



Aug- 



the E+2g perpendicular to the appHed strain (and thus 
experiencing smaller softening), and the E~2g parallel to 
it. This is analogous to the effect of curvature on the 
G peak of carbon nanotubes. The G peak splitting in 
nanotubes is the combined result of electron confinement 
and curvature[H|- Pure curvature splits the graphene E2g 
mode in a component parallel to the tube axis and one 
perpendicular. When the sp^ bonds of graphene are de- 
formed by rolling it in a tube, they lengthen and soften 
in the direction perpendicular to the axis, in order for 
the -Kz electrons to be perpendicular to it. This is pro- 
portional to curvature, so it is minimum parallel to the 
axis, and maximum along the circumference, increasing 
with decreasing diameter [3, 0| . Thus, by curvature only, 



From the data in Fig. 3a we get 7£;2g=1.99; /9E2g=0.99. 
These experimental parameters can now be used to es- 
timate the trends for free-hanging graphene under uni- 
axial strain. Inserting 7£;2g = 1.99, /3£;2g=0-99i i^=0.13 in 
Eq. 4, we get dojc+fde ~-18.6 cm~-^/ %; dioc-lde ~- 
36.4 cm~^/%. Note that the effect of the substrate higher 
Poisson's ratio is to significantly decrease the slope of the 
G+ peak. These results are also in excellent agreement 
with our first-principles calculations (see later). 

We can now use our fitted ^E2g to deduce the expected 
peak variations for graphene under biaxial strain. In 
this case en = eu = £ and, from Eq. 3, AujE^g = 
—2uj%^^^^E2g£^ since the shear deformation term can- 
cels. This means, as expected, that the G peak does 
not split. Also, no difference is expected between free- 
hanging graphene and graphene on a substrate. Thus, 
for biaxial strain: dioolde ~-63 cm~^/ %. 
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To the best of our knowledge, no data exist in literature 
for uniaxial strain on graphite. However, several authors 
applied hydrostatic pressure on graphite 13, 41 , \^ \^ 
finding duJc/dah ~4.4-4.8cm~^/GPa, where c/i is the hy- 
drostatic pressure (stress). The in-plane biaxial strain 
under hydrostatic pressure is e = {Su + Su)(Jh- Since 
for graphite in-plane UiSu + Su) l/1250GPa [sll, 
the data in Refs.fiol. l4ll . 14^ . |43| correspond to an in- 
plane Gruneisen parameter ^E2g ~1.72-1.90, in very good 
agreement with our results. Many groups have consid- 
ered hydrostatic pressure on nanotubes (see.e.g..(33. Issl. 
I43I). Generally it is found dujc/dah ~ 4 — bcm~^ /GPa, 
in good agreement with graphene and graphite. How- 
ever, electron confinement and other effects in nanotubes 
warrant a more detailed comparison of our results on 
graphene with the trends for the individual LO and TO 
G bands in nanotubes, which will be presented elsewhere. 

Several experiments exist for uniaxial strain on 
graphite fibres[4j|. These could be the best approxima- 
tion of uniaxial strain along the graphite plane, since 
their very large diameter compared to single wall nan- 
otubes ensures other possible effects due to electron con- 
finement will be neglig ibleQ. Extensive work on carbon 
fibres of different moduli has shown that the peaks' shift 
is directly related to axial stress rather than strain jiHi. 
Thus, one can assume that in uniaxial experiments the 
applied stress is the known parameter, and the strain ap- 
plied to the atomic bonds can be derived from eu = Sucru , 
where Sii=l/E is the fibre elastic compliance, E the 
Young's modulus and an the applied longitudinal stress. 
Thus, in order to correctly estimate the strain, it is neces- 
sary to know the fibre E, which, in general, is significantly 
lower than the in-plane Young's modulus of graphite |44|. 
Then, if we extend the universal relation between Ra- 
man peak shift and uniaxial stress to graphene, the fol- 
lowing should hold: ^^'^fr^= Eg!^^^^, g"G.^ap,.e„e y^^^^ 
fibres show a uniaxial stress sensitivity of diOGldau ~ 
2 — 3cm~^/GPa In particular, PAN-based carbon 

fibres with "onion skin" morphology (i.e. those most 
similar to large multi-wall nanotubes) have duic / daii=- 
2.3 cm-VGPa Note that, due to disorder, the G 

peak of carbon fibres is very broad and not resolved in 
two subbands. Thus, the fitted G represents the aver- 
age shift of the two subbands. Our average shift, using 
the in-plane graphite Poisson's ratio, as needed in or- 
der to compare with fibres, is dujo/de ~-27cm~^/%. If 
we scale the uniaxial strain sensitivity of PAN fibre s by 
the in plane Young's modulus of graphite~1090GPa [33|, 
this would imply a value of ~-25cm~^/%, in excellent 
agreement with our average value. This also validates 
the assumption that the graphene Young's modulus is 
similar to the in plane Young's modulus of graphite, in 
agreement with recent measurements A notable dis- 
crepancy exists only with Ref.[33| for uniaxial measure- 
ments on fibres. However, their data imply ~2.87 



all graphite literature 0, EH, Si, H SI- We also note 
that our results disagree with recent Raman experiments 



on uniaxial strain in graphene[28|,|32|, which report much 
smaller dui/de, implying much smaller Gruneisen param- 
eters. It is difficult to see how the Gruneisen parameters 
of graphene should be much smaller than those measured 
in-plane for graphite. Moreover, no G peak splitting was 



observed for uniaxial strain again in contrast with 

both our observation and general expectations. 

We now consider the case of the singly degenerate 
phonon modes corresponding to the D and D' peaks. The 
D peak is a breathing mode similar to the TO Aig phonon 
at Kji^l (see Methods). For a pure Aig symmetry and 
small strains, the uniaxial stress shift Aw^j^ is given only 
by the hydrostatic component of the stress: 



-^AigTAig {ett+eu) 



(7) 



On the other hand, the D' phonon is of E symmetry [47| 
and we could expect in principle splitting, and a relation 
similar to Eq. 4. However, experimentally this peak 
is very weak and we cannot resolve any splitting in the 
strain range we have considered. Thus, for small strains, 
we write for both Raman peaks 



AUJ2D: 



2D' 



(8) 



in disagreement with both our measurements and with disorderfl^ [TsT 14 . isl . 2^ . 



Combining our experimental data with Eq. 8 we get 
7£i ~3.55; 7_D' ~1.61. For free-hanging graphene, these 
give duj2D/de ~-83cm-V%; duj2D'/de — 45cm-V% 

In the case of graphene under biaxial strain eu = £44 = 
e and Auj2d,2D' = ^'^^2D:2D''yD;D'£- Thus, using our 
fitted Gruneisen parameters, the expected 2D and 2D' 
variation as a function of biaxial strain are: du)2D/de ^- 
191 cm-V % and duj2D'/de —104 cm'^ / %. 

To the best of our knowledge no data exist for the 2D 
or 2D' peak dependence in grap hite as a function of uni- 
axial strain. However, Ref.[4J| measured duj2Dldcni ~ 
6.4cm~^/GPa for PAN carbon fibres. This scales to 
duj2D/de ~ — 70cTO^^/% in graphene, in agreement 
with our predicted uniaxial trend, when using the in 
plane Possion's ratio of graphite, as needed for compar- 
ison with fibres. For graphite under hydrostatic pres- 
sure Ref.fiill reported duj2D/dah ~12.3 cm~^/GPa, and 
duj2D'ldah ~9 cm~-^/GPa. This corresponds to an in- 
plane biaxial strain e = {Sii+Sit)ah- Vxoml/ {S11+ Sit) ~ 
l/1250GPa Q, we get du2D/de ~ -154cm-V%; 720 = 

2.84; duj2D'/de 113cto-V%; 721)'= 1.74, in broad 

agreement with our predictions for biaxial strain. 

Finally, we note that, in any case, the 2D peak is ex- 
tremely sensitive to strain. With a typical spectrom- 
eter resolution of ~2cm~^, a remarkable sensitivity of 
~0.01% and 0.03% can be achieved for biaxial and uni- 
axial strain, respectively. We also note that a com- 
bined analysis of G and 2D FWHM and shifts should 
allow to dis ting uish between effects of strain, doping or 
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Figure 4: Eigenvectors of G"*" and modes as determined 
by density-functional perturbation theory. Note that these 
are perpendicular to each other, with G~ polarized along the 
strain axis, as expected 



To further understand our findings we perform first- 
principles calculations on free-standing grapliene as de- 
scribed in Methods [131, for small strains up to ^1%, to 
compare with experiments. The effects on electron and 
phonon bands of larger strains will be reported elsewhere. 
Fig. 4 plots the resulting G'^/G~ eigenvectors. These 
are perpendicular to each other, with the G~ eigen- 
vector oriented along the strain direction, as expected. 
For small strains we find dujQ-/de ~-34 cm~^/ % and 
dajQ+fde '^-17 cm~^/%, independent on the strain direc- 
tion, as expected from symmetry. We also get 7E2g = 1.87; 
/3_E2g=0.92, in excellent agreement with our measured pa- 
rameters. Note that, in order to compare the calculated 
trends for G"*" and G~ with our measurements, we need to 
insert the theoretical parameters in Eq. 4 together with 
the substrate Poisson's ratio. This gives dujQ-jde ~- 
30cm~^/ %; dujQ+/de '^-10.3cm~^/%, in excellent agree- 
ment with the fits in Fig. 3a. We also calculate the biaxial 
strain variation for the G peak. We find duoc/de ^-58 
cm~^/%; 7£;2g = 1.8, again in excellent agreement with the 
biaxial values based on our experimental Gruneisen. 

We then calculate the uniaxial and biaxial strain vari- 
ation for the 2D peak. We find du}2D/de ~-60 cm~^% 
for uniaxial, and 81020 1 ~-144 cm~^/ %; for biaxial 
and 7_D ~2.7 for both. These are in excellent agreement 
with the results of hydrostatic pressure experiments on 
graphite, and in broad agreement with our experimen- 
tal data for uniaxial strain (and the consequent biaxial 
predictions), being ~25/% smaller. It is important to 
consider that, while for the Raman active G mode we 
are probing the same centre-zone phonon when measur- 
ing the Raman spectrum on a strained sample, the Ra- 
man D and D' peaks are zone boundary phonons acti- 
vated by double resonance (see Methods). Any change 
in the double Resonance condition during the strain ex- 
periments will vary the actual phonon probed in the 
Raman measurements, as well as inducing a change in 
the phonon frequencies. Thus, the relationship between 
phonon Gruneisen parameters and the variation of the 



Raman peaks with applied strain is in principle more 
complex than the case of the G peak and what implied 
by Eqs 7,8. Indeed, while biaxial strain does not move 
the relative positions of the Dirac cones, uniaxial strain 
changes them [i^. Note that this does not open any 
gap, in contrast with the conclusions of Ref.|28l|. Still, 
it can have a significant infiuence in the double reso- 
nance process. As explained in Methods, while the D' 
is intra-valley, i.e. connecting two points belonging to 
the same cone around K or K', the D peak phonon re- 
quires scattering from the cone around K to that around 
K' 11, 24, [HOl- Thus, its wavevector is determined by 



the relative distance of the Dirac cones and by the laser 
excitation energy. Our experiments are performed for a 
fixed excitation. Then, what we measure in Raman spec- 
troscopy of uniaxially strained graphene is the combina- 
tion of the 2D phonon shift due to strain, and a possible 
additional shift due to the fact that the relative move- 
ment of the Dirac cones changes the phonon wavevec- 
tor. For an asymmetric movement this could lead to 
peak broadening and splitting. Indeed the experimen- 
tal FWHM(2D) significantly increases with strain. In 
the case of the 2D' peak the movement of the relative 
positions of the cones will have no consequence, since 
it is an intra-valley process. However, for both D and 
D', other effects could be given by the renormaHsation 
of Fermi velocity and phonon group velocity with strain. 
Thus, especially for the D peak, our measured has 
to be taken as an upper boundary, and a more general 
expression to evaluate it can be = 



with A'ix'2D encompassing the corrections due to the 
changes in the phonon selected in double resonance, as 
a function of strain. We note that, in the case of biax- 
ial strain, at least the effects due to the relative move- 
ment of the Dirac cones are absent. Thus, Raman ex- 
periments on graphene under biaxial strain would be 
more suited to measure the D mode Gruneisen parame- 
ter, and this explains why our calculations are in excellent 
agreement with the hydrostatic pressure experiments on 
graphite. Thus, given the peculiar nature of electron- 
phonon and electron-electron interactions around the 
K point in graphene[l3, 24, [HH, combined with the 
relative movement of the K, K' points under uniaxial 
strain ji^, and the possible re-normalizations of electron 
and phonon bands, the full theoretical description of the 
2D peak under uniaxial strain needs further investiga- 
tion, and will be reported elsewhere. 

We now consider the polarization dependence of the 
G"*" and G~ intensities, expected due to the nature of the 
phonon eigenvectors and their orientation with respect to 
the strain js^l- The effective photon-phonon interaction 
Hamiltonian for the E2g phonons is |5l| : 

(9) 

Here tx i^y are the cartesian components oi 
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Figure 5: (Color online.) Geometry implied by Eqs. lf9|l- l|lip . 
Tlie circles in the hexagon represent carbon atoms. The x axis 
is chosen to be perpendicular to the C-C bond. The short 
black arrows represent the phonon displacements in the {x, y) 
basis, as assumed in Eq.([9l) (the longitudinal and transverse 
normal modes are given by their linear combinations). The 
strain axis is the blue dashed line. The red arrows represent 
the polarization of incident and detected light. 



the electric field of the incident (scattered) light, and 
Ux,Uy are the phonon displacements in the (a;,y)-basis 
(see Fig. [5] for details). The x axis is chosen perpendic- 
ular to the C-C bond. This Hamiltonian is the only al- 
lowed by the Cgt, symmetry of graphene. In the presence 
of strain the Hamiltonian changes, but the correction will 
be of the order of the strain itself. For a fixed small strain, 
these corrections can be ignored, in first approximation, 
in the calculation of the polarization dependence of the 
G bands. The main effect of strain is to force the phonon 
normal modes to be longitudinal (u/) and transverse {ut) 
with respect to the strain axis, as discussed above, and 
shown in Fig. 4. If we call the angle between the 
strain axis and the x axis, we can write: 

Ux = ui cos ifs -I- Ut sin </Js , Uy ~ ~ui sin + ut cos ips 

(10) 

In our Raman spectrometer, we can excite with lin- 
early polarized light and use an analyzer for the scat- 
tered radiation. This means that the corresponding elec- 
tric field vectors have definite orientations: = 



pin, out 
^0 



cos 



Cq sm( 



where the polarization is measured with respect to the 
strain axis. Substituting these in Eq. ([9]), the matrix 
elements corresponding to emission of longitudinal and 
transverse phonons are proportional to 
3(ps) and sin(6'i„ -I- 9out + 3(p<,), respectively. The intensi- 
ties of the two peaks are given by their squares: 

la- (X cos^{9tn+dout+SVs), Ig+ oc sin^(6'i„+6'o„f+3(y5s) 

(11) 

To test this we do polarization measurements with an 
analyzer for the scattered light aligned with the strain 
direction (flout = 0), and rotating the incident polariza- 
tion with respect to the strain axis in steps of 10°, Fig. 6. 



The data in Fig. 6 are well fitted hy Iq- oc cos 



-56°) 



and Iq+ oc sin^(6'i„ — 56°). According to Eq. (fTT|) . this 
gives ips ~ —18.7°. We thus get the orientation of the 
graphene crystal with respect to the known strain axis. 

The physical origin of the polarization dependence of 
the IG~ peaks can be traced to the microscopic mech- 
anism of Raman scattering. The light interaction with 
graphene phonons is mediated by electrons. As discussed 
in[5l| for unstrained graphene, if one assumes the elec- 
tron spectrum to be isotropic (Dirac), the G peak inten- 
sity vanishes. Thus, the G peak is entirely due to the 
anisotropic terms in the electronic spectrum. In other 
words, in order to contribute to the G peak, electrons 
must "feel" the crystallographic directions. In unstrained 
graphene this has no consequence, since the two vibra- 
tions are degenerate and not resolved. Under strain, the 
two sub-bands correspond to definite orientations of the 
vibrations with respect to the strain axis. It is thus the 
interaction of electrons, which "feel" the crystallographic 
directions, with phonons, entirely determined by the the 
strain direction, that gives the polarization dependence. 

In summary, we probed with Raman spectroscopy the 
optical phonons of graphene as a function of uniaxial 
strain. We find that the doubly degenerate E2g mode 
splits in two components, one polarized along the strain, 
the other, perpendicular. This split of the Raman G peak 
in 2 subbands, G"'" and G~, is analogous to that induced 
by curvature in nanotubes. These subbands red-shift 
with increasing strain, whilst their splitting increases, 
in excellent agreement with first-principles calculations. 
Their relative intensities vary with polarization, allowing 
to probe the sample crystallographic orientation with re- 
spect to the strain. The 2D and 2D' bands downshift, 
but do not split for small strains. Our results can be 
used to quantify the amount of uniaxial or biaxial strain, 
providing a fundamental tool for graphene-based nano- 
electronics and nano/micro electro mechanical systems. 
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METHODS 
Strain and Raman Measurements 

In order to controllably and reproducibly induce strain, 
graphene layers, prepared by micromechanical cleavage 
of graphite, are deposited on two different fiexible sub- 
strates. One is a 720 [mi thick, 23mm long Polyethylene 
terephthalate (PET) film. The other is a 3 mm thick, 10 
cm long, 1 cm wide clear acryHc (Perspex). In both cases 
the large length-to-width ratio is chosen to allow uni- 
form bending and reversibility. Prior to graphene deposi- 
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Figure 6: (Color online.) (Left) Raman spectra and (right) polar plot of the fitted and G~ peaks as a function of the angle 
between the incident light polarization and the strain axis 6i„, measured with an analyzer selecting scattered polarization along 
the strain axis, 6out = 0. The polar data are fitted to Iq- « cos^(Si„ — 56°) and Ig+ oc sin^(0i„ — 56°). 



tion, the substrates are spin coated with SU8 2000.5 (Mi- 
croChem) photoresist [39j of carefully chosen thickness 
(400nm), which is then cross-Hnked. This ensures op- 
timal visible contrast for graphene identification 5^, 5a |. 
To achieve maximum strain, the length of the substrate 
is altered in order to have the fiake at its center. Fig. 
1. Note that the size of the graphene layers is orders 
of magnitude smaller than the substrate length 10'^ 
and ^ 10"* times smaller, respectively). This ensures a 
uniform strain in the section measured by Raman spec- 
troscopy. The first substrate is used in two point bending 
experiments, whilst the second in four point bending. Fig. 
1. Raman spectra are measured with a lOOX objective 
at 514nm excitation with a Renishaw micro- Raman spec- 
trometer, having 1800 grooves/mm grating and spectral 
resolution of ~2cm~^. The polarization of the incident 
light can be controlled by a Fresnel rhomb, while an ana- 
lyzer can be placed before the grating. The power on the 
samples is well below 2mW, so that no shift, nor change in 
width of the Raman peaks is observed for a fixed strain, 
thus ensuring no damage, nor heating. A cycle of load- 
ing, unloading and loading is followed to ensure repro- 
ducibility for both experiments. A total of 80 Raman 
spectra are measured for an average strain increment of 



0.05%. The maximum strain appHed to the sample is less 
than ~1.2%. In the two point measurements, the spec- 
tra do not change until a nominal strain of ~0.55% is 
applied to the substrate. Afterwards they evolve linearly 
with strain. Thus, we assume this point as the reference 
zero strain for the sample. In the four point, the spec- 
tra evolve linearly from zero strain. The two set of data 
are fully overlapping, further confirming the strain mea- 
surements. The data are fully reproducible over three 
strain cycles between maximum and minimum, as shown 
in Fig. 3. Only when suddenly applying large strains or 
large strain increments we observe sample sHppage, indi- 
cated by an upshift, or smaller downshift or no shift at all 
of the Raman parameters. Indeed, for samples suddenly 
bent to large strain values of a few % we often observe 
no change in the Raman peaks, indicating a general loss 
of contact between the graphene and the substrate. It is 
thus extremely important to apply the strain in the most 
controlled way in order to ensure reproducibility and no 
slippage. A further set of 36 measurements is done for a 
fixed value of strain, by rotating the incident polarization 
in 10° steps with respect to the strain axis, and analyzing 
the scattered fight in the plane parallel to the strain axis. 
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Origin of the Raman peaks 



All carbons show common features in their Raman 
spectra in the 800-2000 cm~^ region, the so-called G 
and D peaks, which He at around 1580 and 1350 cm~^ 
respectively [5J|. The G peak corresponds to the dou- 
bly degenerate E2g phonon at the Brillouin zone center. 
The D peak is due to the breathing modes of sp^ rings 
and requires a defect for its activation 5^, H^l- It comes 
from TO phonons around the K point of the Brillouin 
zone [551, [501, is active by double resonance 



and is strongly dispersive with excitation energy due 
to a Kohn Anomaly at Kji^l- The activation process 
for the D peak is an inter- valley process as follows: i) 
a laser induced excitation of an electron/hole pair; ii) 
electron-phonon scattering with an exchanged momen- 
tum q K; iii) defect scattering; iv) electron/hole re- 
combination. The D peak intensity is not related to the 
number of graphene layers, but only to the amount of 



disorder 
cess i.e. 



55 



56]. DR can also happen as intra- valley pro- 
connecting two points belonging to the same 
cone around K (or K'). This gives rise to the so-called 
D'peak, which can be seen around 1620 cm~^ in defected 
graphite [5§|. The 2D peak is the second order of the 
D peak. This is a single peak in monolayer graphene, 
whereas it splits in four bands in bilayer graphene, re- 
flecting the evolution of the band structure [11|. The 2D' 
peak is the second order of the D' peak. Since 2D and 2D' 
peaks originate from a process where momentum conser- 
vation is obtained by the participation of two phonons 
with opposite wavevectors (q and — q), they do not re- 
quire the presence of defects for their activation, and are 
thus always present. Indeed, high quality graphene shows 
the G, 2D and 2D' peaks, but not D and D'[iH. 



Density Functional Calculations 



We use density-functional theory (DFT ) an d density- 
functional perturbation theory (DFPT) [59(1 as imple- 
mented in the PWSCF package of the Quantum- 
ESPRESSO distribution (H^l, within the local-density 
approximation [gH, with norm-conserving pseudopoten- 
tials [g^I and a plane-wave expansion up to 55 Ry 
cut-off. The Brillouin-zone is sampled on a 42x42x1 
Monkhorst-Pack mesh for graphite and graphene, with a 
cold smearingje^l in the electronic occupations of 0.02 Ry. 
We use the equilibrium lattice parameter a = 2.43 A and 
an interlayer spacing of 15 A. We apply the strain in 
different directions. For each direction and strain we de- 
termine the structure with the lowest total energy, by 
varying the size of the unit cell in the direction perpen- 
dicular to the strain. Our calculated values at zero strain 
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